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Abstract 

We describe a new class of scattering matrices for quantum graphs in which back- 
scattering is prohibited. We discuss some properties of quantum graphs with these 
scattering matrices and explain the advantages and interest in their study. We also 
provide two methods to build the vertex scattering matrices needed for their construc- 
tion. 

1 Introduction 

It has been proposed that quantised versions of metric graphs could be used to investigate 
the origin of spectral correlations in semi-classical quantum systems pj. This programme 
has borne a considerable amount of fruit in the past decade. The culmination has been 
the ground-breaking work [2] providing a mechanism to understand the fidelity of spectral 
correlations to random matrix theory in a large class of so-called quantum graphs. It is now 
clear that quantum graphs can be considered at the forefront of attempts to understand 
the universal spectral correlations present in quantised chaotic systems [3l H] . For a review 
with ample references we refer the reader to [5]. 

One interesting new development has been the discovery that the techniques of [T] 
can be used to develop trace formula and associated spectral zeta functions for discrete 
Laplacians on graphs [6j. This leads to the possibility of investigating deep connections 
between quantum graphs and combinatorial spectral graph theory, hitherto unexploited. 
In the present article we expand on some ideas introduced in [6]. 

We consider finite, connected graphs without loops or multiple bonds, and which have 
B bonds and V vertices. The topology of the graph is specified in terms of the V x V 
connectivity matrix C, whose p, q element takes the value 1(0) if the vertices p and q are 
connected (not connected). We denote by vj the valency (degree) of the vertex j. A graph 
for which all the vj are equal to v, is called a f-regular graph. The metric properties are 
provided by the bond lengths Lf, , 1 < b < B. 

Spectral problems on quantum graphs are often written in terms of a 2B x 2B unitary 
quantum evolution matrix U defined on the vector space spanned by the directed bonds [6] . 
U itself can be defined in terms of the vj x vj unitary scattering matrices a^^\ j = 1 . . .V 



associated with each vertex of the graph. 



where (pj) is a directed bond with initial vertex p and terminal vertex j. The spectrum 
under consideration is the set of values {kn}nGZ of the spectral parameter k for which 

det[/2B - Uik)] = 0. (1.2) 

(Here, and elsewhere denotes the identity matrix of dimension n.) A vertex scattering 
matrix a may be derived from boundary conditions of a self-adjoint operator as in [T], or 
specified a priori in order to provide a wider class of examples [TJ [8] , which is the approach 
we adopt here. 

The classical analogue of a quantum graph is a Markovian process on the directed 
bonds [T] with matrix M of transition probabilities related to U by 

T\yr \TT 1 2 

^^J-{pj){mq) — \^{p'j){mq)\ ■ 

Thus M is a doubly stochastic matrix, and hence it has dominant eigenvalue 1. It can be 
interpreted as the classical operator which induces a discrete time Markovian evolution on 
the graph. The spectral gap A is the gap between the leading eigenvalue 1 and the second- 
largest eigenvalue of M, in absolute value. It is conjectured that the rate of convergence of 
spectral statistics of quantum graphs to random matrix theory is governed by the rate of 
ergodicity in this Markov evolution. More precisely, for a sequence of graphs with B ^ oo 
the eigenvalues of the quantum graph are given by random matrix theory if BA ^ cxd |8]. 

We propose to study a new class of vertex scattering matrices defined in the following 
way. 

Definition 1.1. A v x v unitary matrix a is equi-transmitting if (Tpp = for all p, and 
the off-diagonal elements have equal amplitudes; \(Tpq\ = {v — 1)~2 for p ^ q, where v is 
the valency of the vertex. 



Our definition of an equi-transmitting matrix involves two properties which have impli- 
cations for the classical analogue of the quantum graph. The first ensures that a classical 
particle on a directed bond has probability zero to be back-scattered to the reversal of 
that bond in the next time step. The second property produces democratic transmission 
probabilities M(pj)(^g) = {vj - l)~^5mj, for p / g 

Other common choices for vertex scattering matrices are the Neumann matrix (t'^I; 

4^1 = ^ - S„, (1.3) 

and Fourier transform matrices o"!^' introduced as vertex scattering matrices in [8], 

= -Le^-'P"/- . (1.4) 

Fourier transform scattering matrices allow back-scattering with equal transmission and 
reflection probabilities, M{j,j){mq) = v~^5mj- The Neumann scattering matrices prefer 
back-scattering transitions over all other transitions put together. 
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The transition probabilities induced by the three classes of scattering matrices at a 
vertex with a valency v, can be summarised by 

|o-iip = r and \aij\'^ ^ - 1 ^^^^'^^^ (^•^) 

where r = (2/t>— 1)^ (Neumann), r = 1/v (Fourier transform) and r = (equi-transmitting) . 
The aims of the present manuscript can be summarised as follows: 

1. To explain the motivations for studying graphs with equi-transmitting vertex scat- 
tering matrices. 

2. To prove that equi-transmitting vertex scattering matrices exist for arbitrarily large 
V, and to provide practical methods to construct them. 

3. To investigate their spectral properties on both the classical (gap estimates) and the 
quantum (spectral statistics) levels. 

The next section will explain the motivation for introducing the equi-transmitting 
scattering matrices. While equi-transmitting matrices will be shown to have desirable 
qualities, an elementary calculation shows it is trivial to construct 2x2 examples, but 
no 3 X 3 equi-transmitting matrix exists. In section [3] we address the question of their 
existence in other dimensions, and provide several infinite sets of examples. Whether 
equi-transmitting matrices exist in all dimensions greater than three is an interesting 
open question. In section U] we compare spectra of the matrix M for graphs quantised 
with different scattering matrices. We find that in some fairly general situations the 
spectral gap in equi-transmitting quantum graphs is larger than that obtained with other 
scattering matrices. In section[5]we present the results of numerical simulations using equi- 
transmitting quantum graphs that show that random matrix statistics are reproduced. 



2 Motivations 

In this section we shall try to explain the motivations for the introduction of the equi- 
transmitting vertex matrices. 

A fundamental building block in the spectral theory for quantum graphs and a key 
tool in understanding spectral correlations [5], is the trace formula [HlH]. It relates the 
quantum spectrum of a graph, {kn}nez, to the length spectrum of its periodic orbits, 
expressed as an identity of distributions: 

^ e*'^"" = ^5iu) + - ^p) + ^P'^(" + ^p)) • (2-6) 

nGZ P-o. P 

The first term on the right-hand side of (j2.6p is the Weyl term. ^ denotes twice the 
total length of the graph, and plays the role of the volume. The more interesting second 
term is a sum which goes over classes of closed itineraries on the graph equivalent up to 
cyclic permutations of the edges. We call these periodic orbits of the graph. The metric 
length of an orbit is denoted ip, and Vp is the number of times the orbit is a repetition 
of a shorter one. The amplitude factor Ap is the product of all elements of the vertex 
scattering matrices encountered as the orbit passes from one bond to the next. 
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Since equi-transmitting matrices forbid back-scattering transitions, any orbit in which a 
traversal of a bond is fohowed immediately by the traversal of its reverse is eradicated from 
the sum. This significantly reduces the number of orbits which need to be considered — in 
fact the orbits that remain are exactly the closed geodesies considered in combinatorial 
graph theory which will be discussed below. For a ?;-regular graph the asymptotic number 
of orbits of period n is reduced from w"/n to {v — l)"/n. 

We anticipate that using equi-transmitting scattering matrices will significantly sim- 
plify periodic orbit theories in quantum graphs, and perhaps lead to new interesting prob- 
lems and breakthroughs. The following will serve as an example: A theorem of Gutkin and 
Smilansky [lOj guarantees that one can "hear" the shape of a graph if the bond lengths 
are rationally independent and if the vertex scattering matrices are "properly connecting" . 
The equi-transmitting scattering matrices do not belong to the latter class, and the ques- 
tion if one can "hear" equi-transmitting graphs is open. The proof of the above mentioned 
theorem is based on the trace formula, and relies heavily on the special properties of the 
2-periodic orbits which are absent from the trace formula of equi-transmitting graphs. 
Thus, a completely novel approach has to be developed. 

The study of quantum graphs with equi-transmitting vertex scattering matrices leads 
to connections with objects that have been extensively studied in combinatorial graph 
theory, such as the Ihara-Selberg zeta function [11] and its generalisations. One way to 
see this connection is by studying the spectrum of the classical evolution operation M, 
which is defined as the zero set of the secular function 



This equation can be used to derive a classical (Rouelle like) trace formula which is based 
on the periodic orbits on the graph weighted by the products of scattering probabilities 
along the orbit. When dealing with equi-transmitting vertex scattering matrices, the 
trace formula (and the corresponding zeta function) includes only periodic orbits without 
back-scatter. These are the orbits which appear in the Ihara-Selberg zeta function and 
its extensions [12^ [T3l [14} I15j . Moreover, the zeta functions are related by simple trans- 
formations, so that an interesting correspondence between the two seemingly unrelated 
problems can be established. We have mentioned previously the classical spectral gap 
A, and its conjectured influence on the spectral statistics in the corresponding quantum 
graph. Establishing bounds on A arises in the study of the Ihara-Selberg zeta function, 
which makes an intriguing link between spectral statistics and number theory. (See [16] 
for a review of the connections between arithmetical zeta functions and spectral statistics.) 
We shall make use of these connections here (see section H]) . 

We are not able to prove that graphs with equi-transmitting vertex scattering matrices 
display spectral statistics which reproduce the predictions of random matrix theory. We 
present, however, quite convincing numerical evidence showing that both their spectral 
repulsion and spectral rigidity adhere to the predictions of the canonical random matrix 
ensembles. Whether this observation is valid and can be rigorously formulated is an open 
problem awaiting future research. 



Zm{h) :=det[/2B -/uM] . 
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3 Existence and construction 



In this section we shall show that the set of equi-transmitting matrices is not empty or 
trivial. It is easy to construct equi-transmitting matrices in dimension 2. It is equally 
easy to show that equi-transmitting matrices do not exist in dimension 3. We are not 
able to provide a list of the dimensions for which equi-transmitting matrices exist. We 
can, however, show that this list is infinite. We do it by constructing examples of equi- 
transmitting matrices using skew-Hadamard matrices [n\ \TE[ [T9] , and Dirichlet characters 

m- 

3.1 Construction of equi-transmitting matrices using Hadamard 
matrices 

Definition 3.1. A Hadamard matrix H is a matrix whose entries are ±1, and whose 
columns are orthogonal. A skew-Hadamard matrix is a Hadamard matrix satisfying the 
additional condition 

H + H^ = 2/^, 



Proposition 3.2. Let H be a v x v skew-Hadamard matrix. Then 



a = ^L={H-I,) (3.8) 



is an equi-transmitting matrix. 



Proof. All entries of a are ib(u — 1) except for the zero entries along the diagonal. So 
we need only check unitarity. But this is clear, since 

CTCT^ = {HH^ -{H + H^) + h) 

= Iv (3.9) 

since necessarily HH^ = vly. □ 

Hadamard matrices have been conjectured to exist in dimensions 1,2 and all multi- 
ples of 4. This conjecture appears to date back to [21]. Various constructions of skew- 
Hadamard matrices are known [n\ [T8l I19j . Currently there are constructions of skew- 
Hadamard matrices for all dimensions which are a multiple of 4 up to and including 184 
[221 (23], plus other infinite sets of dimensions. 



3.2 Construction of equi-transmitting matrices using Dirichlet 
characters 

This method provides equi-transmitting matrices of dimensions P + 1 where P is prime. 
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Proposition 3.3. Let P he an odd prime and let x a non-trivial Dirichlet character 
modulo P. Let C he the P x P matrix defined hy Cji = ~ j)- Then 



/ 1 ••• 1 \ 



a 



c 



\ 1 



(3.10) 



is equi-transmitting. 



We first need an auxiliary lemma. 
Lemma 3.4. Let P be an odd prime, x non-trivial Dirichlet character, and j £ Z. Then 

P-i 



^ X{n- j)x{n) 



n=0 



P-l ifj = modP, 
— 1 otherwise. 



(3.11) 



Proof. If j = modulo P then 

P-i P-i 
5;ix(n)P = j;i = P-l. 

n=0 n=l 

Otherwise, writing for the multiplicative inverse of n in the finite field 

p-i p-i 



(3.12) 



X{n - j)x{n) = '^x{jn- j)x{n) 

n=0 n=l 
P-l 



-1 



^X{{n-j)n ^] 



n=l 
P-l 



jn 



n=l 



(3.13) 



Now since n ^ runs over all invertible elements in the field, the only argument which does 
not appear in the sum is 1. So 

p-l p-l 

X{n - j)xin) = X{m) - x(l) 
n=0 m=0 

= -1, (3.14) 

Using orthogonality of Dirichlet characters [20] and the known value x(l) = 1- D 

Proof of proposition 's^ We first observe that |x("')l = 1 unless n = whence x(0) = Oi 
so a has the required form. It is easy to see that unitarity will follow once we prove that 



CC^ 



f P-l -1 
-1 P-l 

V -1 -1 



-1 \ 

-1 

p-l J 



(3.15) 
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For this, note that the inner product of the j*^ and i^^ columns of C can be written 

p 



m=l 



P 

E 

m=l 
P-1 

E 

n=0 



X(j - m)x{^ - m) 



X{n + j - i)x{n), 



(3.16) 



via a change of index of summation. Now (|3.15|) follows by using the lemma to evaluate 
the sum (f3l^ . □ 

Corollary 3.5. Let P be a prime congruent to 1 modulo 4. Then there exists a symmetric 
equi-transmitting matrix of dimension P + 1. 

Proof. We use the construction in proposition 13. 31 with the Legendre symbol as the Dirich- 
let character xi^) — (t) '^^^re 

if n = mod P, 

1 if n is a square modulo P, (3.17) 
— 1 if n is not a square modulo P. 



To show that a is symmetric it suffices to show that the circulant matrix C is symmetric. 
To see this, note that 



xiP 



m 



m{ — l) 
P 



-1 



p) 



X{m) 



since (-1/P) = 1 if P = 1 modulo 4 (Euler). □ 

By Dirichlet's theorem, corollary 13.51 provides infinitely many examples of symmetric 
equi-transmitting matrices. 

The constructions in propositions 13.21 and 13.31 give many examples of equi-transmitting 
matrices of even dimensions. We can construct an example of an equi-transmitting matrix 
in dimension 5, 



1 

a = - 
2 





1 


1 


1 


1 


\ 


1 





1 


UJ 






1 


1 











1 


U! 







1 




V 1 






1 





/ 



where oj = e^'^'/^. 



(3.18) 



However, apart from this example we do not have any examples in odd dimensions. These 
equi-transmitting matrices do not appear to have been studied in the literature. An 
examination of (|3.18p might lead one to suspect that such matrices can be constructed for 
dimension v by using entries that are (f — 2)*^^ roots of unity. However we have exhaustively 
checked for the case v = 7 and shown this to be false. 



4 Properties of graphs with equi-transmitting scattering 
matrices 

Now that we have demonstrated the existence of equi-transmitting matrices, we can study 
the quantum and classical evolutions which they induce. In particular, we shall use some 
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results from combinatorial graph theory to demonstrate the advantages gained by studying 
graphs with equi-transmitting vertex scattering matrices. 

An important tool in the preceding discussion is the spectrum of the connectivity 
matrix C defined in the introduction. Since C is symmetric its eigenvalues /ij are real, 
and we order them: ^v-i ^ • • • ^ /^i ^ A^o- 

We shall consider u-regular graphs and to avoid trivial cases we will assume throughout 
that f > 3. A connected f-regular graph has the property that its connectivity matrix 
has largest eigenvalue fiQ = v and it is simple. A u-regular graph is called Ramanujan if 
all other eigenvalues of C are contained in the interval [— 2-v/t7^^, 2\/t7^^]. Ramanujan 
graphs are of interest in computer science and communication network theory since they 
are sparse yet highly connected [25] . 

We shall discuss the spectrum of the classical evolution operator M on a w-regular 
graph constructed by using vertex scattering matrices of the types listed in (jl.Sp . We 
shall relate them to the spectrum of the connectivity matrix C by the following theorem. 

Theorem 4.1. Let M be the doubly stochastic transition probabilities matrix associated 
to a quantum v-regular graph with unitary vertex scattering matrices a which satisfy (jl.5|) 
for some r > 0. Let the eigenvalues of C be fiQ ^ fii ^ ■ ■ ■ ^ l^v-i- Let 



(1 - r)fij + ^{1 - rYiJ?- - 4(1 - rv){v - 1) 



2{v - 1) 



and 



(1 - r)yij - ^{l-rYpL]-A{l-rv){v-l) 



2{v-l) 



The spectrum of M consists of the points 

|1 — rv\ |1 — rv\ 



Uq,... ,UV-1,Uo, . . .,UV-1, 



V — 1 V — 1 

{v - 2)V 



where the last two points are listed with multiplicity 

z 

Proof. Let W be the matrix representing Hashimoto's bond (edge) adjacency operator 
[Ml IS] defined as: 

^{pj)(m<?) = Smji^ - Spq). (4.19) 

It has entries equal to 1 only when two directed bonds follow each other at a common 
vertex j, but excluding back-scattering. The form of the vertex scattering matrices cr of 
interest here is provided by (jl.Sp . It implies that for w-regular graphs, 

where r = (2/f — 1)^ (Neumann), r = 1/v (Fourier transform) and r = (equi-transmitting). 
The characteristic polynomial of M, det[ul2B — M] is related to a graph theoretic zeta 
function developed in [27J. Bartholi's theorem \27\ I28j implies an equivalent form for the 
characteristic polynomial in terms of the matrix C, 



(I - r )2\ (''"2)^2 
det[uL2B - M] = ( - _Y^2 J 



+ Lv -Cu 

V — 1 J V — 1 
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It follows that eigenvalues of M are solutions to 



(1 - rv)^ ^ 
(v-l)^ ' 

(with multiplicity (f — 2)V/2) and 

2 1 — r 1 — rv 

u -fj,ju H — = 0, 

V — 1 V — 1 

for J = 0, . . . , F - 1. n 

For regular graphs with equi-transmitting scattering matrices, theorem 14. II shows that 
the eigenvalues of M are (up to a scaling) at the positions of the poles of the Ihara-Selberg 
zeta function [llj of the graph, as was noted in [6]. The proof in this case follows from 
Bass' identity [29j of which Bartholi's theorem is a generalisation. 

Theorem 14.11 will enable us to compare equi-transmitting scattering matrices with 
others of type (jl.Sp (see Theorem 14.31 below). It also has a number of other consequences 
which may be of independent interest. 

Corollary 4.2. • Ifr^l/v then all eigenvalues of M are real. 

• If r = \/v (e.g. Fourier transform scattering matrix) then the eigenvalues of M are 

; • • • 5 : 

V V 

and the has multiplicity {v — 1)V . 



,0 



2 



If r < 1/v then Uj and Uj are real iff ^ ^/ {1 — rv){v — 1). 

2 

If r < 1/v and \fij\ < y (1 — rv){v — 1) then 



Ui 



rv 



1 



In particular, for a u-regular graph which is Ramanujan, and has equi-transmitting 
scattering matrices, all but the eigenvalue 1 lie in a disc of radius {v — 1)"^/^ about the 
origin. 

Theorem 4.3. Consider a v-regular graph. Let the spectral gap for the quantum graph 
with equi-transmitting scattering matrices be Act and denote by the spectral gap for 
the same graph with scattering matrices a satisfying ()1.5p for some r > 0. Then there 
exists e > such that if C has an eigenvalue in either of the intervals (2\/v — 1 — e, v) or 
[—V, —2^v — 1 + e) then Aet > A^. 

We give the proof of theorem 14.31 at the end of this section. Theorem 14.31 demonstrates 
that in some fairly general situations the spectral gap arising in quantum graphs with 
equi-transmitting scattering matrices is larger than the spectral gap with other kinds of 
scattering matrices. Notice that if one considers a sequence of f-regular graphs with 
y — > oo the Alon-Boppana bound Theorem 1.3.1] states that 



liminf ui > 1\/v — 1. 
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In other words, theorem 14.31 will apply eventually. 

In the interests of full disclosure we point out that for some other families of graphs equi- 
transmitting matrices will not necessarily lead to a larger (although still large) spectral 
gap, e.g. fully-connected graphs. 

Before giving the proof, we provide one more lemma. 
Lemma 4.4. Let ^ < v. Then define 



, . _ (1 - r)fi + ^{l - rfn"^ + A{vr - l){v ^ 
^^^''^ 2{v - 1) 

Then f^{r) is real and strictly increasing on 1 ^ — - — -(v — \J — /i^) ^ r ^ 1. 

Proof. It is convenient to define a new variable X = v{l — r)/{v — l), so that becomes, 
after some manipulation, 

and is real if ^ X ^ — — ^^^^ ^ - . Differentiating, 

dU _ 1 2v^-i?X \ 



< — 



1 / {2-X)v 



2 



2V \^ y/w2(X-2)2 J 

= i-(/.-.)<0 

since ^/{X - 2)2 = 2 - X, as ^ X < 2. As ^ < it follows that /^(r) is strictly 
increasing. □ 

Proof of theorem \4.3[ We assume ^ l/iy-il. If this is not the case then the argument 
below holds mutatis mutandis replacing fii by nv-i- 

If the graph is Ramanujan, and if r < 1/v, then choose e so that 

^—^Vv^<2V^irn:-e. (4.20) 
1 — r 

Since f > 3 if r > such an e can always be found. 

2 — rv 2 

As fii > \/v — 1 > Y^(l — rv)(v — 1) by corollary 14.21 ui = fnAr) is real 

1 — r 1 — r 
and Aj. = 1 — ui, where 



_ (1 - r)^i + V(l - r)2^2 _ _ ^^)(^ _ 1) 1 
^ 2(t;-l) V^T^' ^ ■ ^ 

For Ramanujan graphs, Aet = 1 1= > A,.. 

Vf — 1 
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If the graph is Ramanujan and r ^ 1/v then choose e to satisfy 

<2^/v-l-e. (4.22) 

Vf - 1 

By lemma 1331 

u,^f,,{l/v) = ^>^L=, (4.23) 

and again we have Act > ^r- 

If the graph is not Ramanujan, then ui is still real since 

2 

III > 2Vf -l> -\/ (1 - rv){v - 1) ifr<l/u, 

so by lemma 133] 

1 - A, = f^, (r) > f^, (0) = 1 - Aet. (4.24) 

□ 

5 Numerical simulations 

The results of the previous section show that well-connected quantum graphs with equi- 
transmitting scattering matrices can have large spectral gaps. According to the Tanner 
conjecture [8] the spectral statistics should converge to the statistics of ensembles of ran- 
dom matrix theory. To illustrate this we present the results of some numerical calculations 
of the nearest-neighbour spacing density P{s) for points in the spectrum of a quantum 
graph, and the variance V{L) of the number of points in an interval of length L. 

To calculate these spectral statistics we did not solve (|1.2p directly. Rather we took 
an approach which is known to be equivalent. We replaced the phases fe-^^(pj) in (ll-ip 
by random phases in the interval [0, 27r] and diagonalised the resulting matrices. The 
statistics of the re-scaled eigen-phases approach those of the spectrum defined by (|1.2p 
in the limit as i? ^ oo if the bond lengths are not rationally related, and drawn from a 
narrowing interval as — > oo O [30] . 




Figure 1: The graphs used in the numerical simulations: a) A complete graph with 13 
vertices, b) a 5-regular graph with 20 vertices. 
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In figure [2] we plot the nearest neighbour distribution and number variance for the 
5-regular graph in figure [TJ)) with the equi-transmitting vertex scattering matrix given 
exphcitly in ()3.18|) . Since a is symmetric we expect the spectral statistics to approach 
those of the Gaussian orthogonal ensemble (GOE) of random matrices. We also plot in 
figure [2] the corresponding limiting curves for Gaussian orthogonal ensemble and Gaussian 
unitary ensemble (GUE) [31j. We see, as expected, agreement to the Gaussian orthogonal 
ensemble curves, even for such a relatively small graph. 

In figure[3]we present the corresponding numerics for the complete graph on 13 vertices 
(figure la) with an anti-symmetric equi-transmitting scattering matrix at each vertex. In 
this case we expect convergence to the Gaussian unitary ensemble statistics, and this is 
also clearly demonstrated in the figure. 



Pis 




Figure 2: Spectral statistics for a quantum regular graph with symmetric scattering ma- 
trices. On the left is a plot of the nearest-neighbour density; on the right is the number 
variance. 



Given the utility of these equi-transmitting matrices, we believe the study of their 
existence in dimensions for which we do not currently have examples is interesting and 
merits further investigation. 
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Figure 3: Spectral statistics for the quantum complete graph with anti-symmetric scat- 
tering matrices. On the left is a plot of the nearest-neighbour density; on the right is the 
number variance. 
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